skip to main content


Search for: All records

Creators/Authors contains: "Johnson, James W."

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. ABSTRACT

    The scaling of the specific Type Ia supernova (SN Ia) rate with host galaxy stellar mass $\dot{\text{N}}_\text{Ia} / \text{M}_\star \sim \text{M}_\star ^{-0.3}$ as measured in ASAS-SN and DES strongly suggests that the number of SNe Ia produced by a stellar population depends inversely on its metallicity. We estimate the strength of the required metallicity dependence by combining the average star formation histories (SFHs) of galaxies as a function of their stellar mass with the mass–metallicity relation (MZR) for galaxies and common parametrizations for the SN Ia delay-time distribution. The differences in SFHs can account for only ∼30 per cent of the increase in the specific SN Ia rate between stellar masses of M⋆ = 1010 and 107.2 M⊙. We find that an additional metallicity dependence of approximately ∼Z−0.5 is required to explain the observed scaling. This scaling matches the metallicity dependence of the close binary fraction observed in APOGEE, suggesting that the enhanced SN Ia rate in low-mass galaxies can be explained by a combination of their more extended SFHs and a higher binary fraction due to their lower metallicities. Due to the shape of the MZR, only galaxies below M⋆ ≈ 3 × 109 M⊙ are significantly affected by the metallicity-dependent SN Ia rates. The $\dot{\text{N}}_\text{Ia} / \text{M}_\star \sim \text{M}_\star ^{-0.3}$ scaling becomes shallower with increasing redshift, dropping by factor of ∼2 at 107.2 M⊙ between z = 0 and 1 with our ∼Z−0.5 scaling. With metallicity-independent rates, this decrease is a factor of ∼3. We discuss the implications of metallicity-dependent SN Ia rates for one-zone models of galactic chemical evolution.

     
    more » « less
  2. ABSTRACT

    We model the stellar abundances and ages of two disrupted dwarf galaxies in the Milky Way stellar halo: Gaia-Sausage Enceladus (GSE) and Wukong/LMS-1. Using a statistically robust likelihood function, we fit one-zone models of galactic chemical evolution with exponential infall histories to both systems, deriving e-folding time-scales of τin = 1.01 ± 0.13 Gyr for GSE and $\tau _\text{in} = 3.08^{+3.19}_{-1.16}$ Gyr for Wukong/LMS-1. GSE formed stars for $\tau _\text{tot} = 5.40^{+0.32}_{-0.31}$ Gyr, sustaining star formation for ∼1.5–2 Gyr after its first infall into the Milky Way ∼10 Gyr ago. Our fit suggests that star formation lasted for $\tau _\text{tot} = 3.36^{+0.55}_{-0.47}$ Gyr in Wukong/LMS-1, though our sample does not contain any age measurements. The differences in evolutionary parameters between the two are qualitatively consistent with trends with stellar mass M⋆ predicted by simulations and semi-analytic models of galaxy formation. Our inferred values of the outflow mass-loading factor reasonably match $\eta \propto M_\star ^{-1/3}$ as predicted by galactic wind models. Our fitting method is based only on Poisson sampling from an evolutionary track and requires no binning of the data. We demonstrate its accuracy by testing against mock data, showing that it accurately recovers the input model across a broad range of sample sizes (20 ≤ N ≤ 2000) and measurement uncertainties (0.01 ≤ σ[α/Fe], σ[Fe/H] ≤ 0.5; $0.02 \le \sigma _{\log _{10}(\text{age})} \le 1$). Due to the generic nature of our derivation, this likelihood function should be applicable to one-zone models of any parametrization and easily extensible to other astrophysical models which predict tracks in some observed space.

     
    more » « less
  3. ABSTRACT

    We derive empirical constraints on the nucleosynthetic yields of nitrogen by incorporating N enrichment into our previously developed and empirically tuned multizone galactic chemical evolution model. We adopt a metallicity-independent (‘primary’) N yield from massive stars and a metallicity-dependent (‘secondary’) N yield from AGB stars. In our model, galactic radial zones do not evolve along the observed [N/O]–[O/H] relation, but first increase in [O/H] at roughly constant [N/O], then move upward in [N/O] via secondary N production. By t ≈ 5 Gyr, the model approaches an equilibrium [N/O]–[O/H] relation, which traces the radial oxygen gradient. Reproducing the [N/O]–[O/H] trend observed in extragalactic systems constrains the ratio of IMF-averaged N yields to the IMF-averaged O yield of core-collapse supernovae. We find good agreement if we adopt $y_\text{N}^\text{CC}/y_\text{O}^\text{CC}=0.024$ and $y_\text{N}^\text{AGB}/y_\text{O}^\text{CC} = 0.062(Z/Z_\odot)$. For the theoretical AGB yields we consider, simple stellar populations release half their N after only ∼250 Myr. Our model reproduces the [N/O]–[O/H] relation found for Milky Way stars in the APOGEE survey, and it reproduces (though imperfectly) the trends of stellar [N/O] with age and [O/Fe]. The metallicity-dependent yield plays the dominant role in shaping the gas-phase [N/O]–[O/H] relation, but the AGB time-delay is required to match the stellar age and [O/Fe] trends. If we add ∼40 per cent oscillations to the star formation rate, the model reproduces the scatter in the gas phase [N/O]–[O/H] relation observed in external galaxies by MaNGA. We discuss implications of our results for theoretical models of N production by massive stars and AGB stars.

     
    more » « less
  4. Abstract

    We measure abundances of 12 elements (Na, Mg, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni) in a sample of 86 metal-poor (−2 ≲ [Fe/H] ≲ −1) subgiant stars in the solar neighborhood. Abundances are derived from high-resolution spectra taken with the Potsdam Echelle Polarimetric and Spectroscopic Instrument on the Large Binocular Telescope, modeled using iSpec and MOOG. By carefully quantifying the impact of photon-noise (<0.05 dex for all elements), we robustly measure theintrinsicscatter of abundance ratios. At fixed [Fe/H], the rms intrinsic scatter in [X/Fe] ranges from 0.04 (Cr) to 0.16 dex (Na), with a median of 0.08 dex. Scatter in [X/Mg] is similar, and accounting for [α/Fe] only reduces the overall scatter moderately. We consider several possible origins of the intrinsic scatter with particular attention to fluctuations in the relative enrichment by core-collapse supernovae (CCSN) and Type Ia supernovae and stochastic sampling of the CCSN progenitor mass distribution. The stochastic sampling scenario provides a good quantitative explanation of our data if the effective number of CCSN contributing to the enrichment of a typical sample star isN∼ 50. At the median metallicity of our sample, this interpretation implies that the CCSN ejecta are mixed over a gas mass ∼6 × 104Mbefore forming stars. The scatter of elemental abundance ratios is a powerful diagnostic test for simulations of star formation, feedback, and gas mixing in the early phases of the Galaxy.

     
    more » « less
  5. Abstract We investigate the [X/Mg] abundances of 16 elements for 82,910 Galactic disk stars from GALAH+ DR3. We fit the median trends of low-Ia and high-Ia populations with a two-process model, which describes stellar abundances in terms of a prompt core-collapse and delayed Type-Ia supernova component. For each sample star, we fit the amplitudes of these two components and compute the residual Δ[X/H] abundances from this two-parameter fit. We find rms residuals ≲0.07 dex for well-measured elements and correlated residuals among some elements (such as Ba, Y, and Zn) that indicate common enrichment sources. From a detailed investigation of stars with large residuals, we infer that roughly 40% of the large deviations are physical and 60% are caused by problematic data such as unflagged binarity, poor wavelength solutions, and poor telluric subtraction. As one example of a population with distinctive abundance patterns, we identify 15 stars that have 0.3–0.6 dex enhancements of Na but normal abundances of other elements from O to Ni and positive average residuals of Cu, Zn, Y, and Ba. We measure the median elemental residuals of 14 open clusters, finding systematic ∼0.1–0.4 dex enhancements of O, Ca, K, Y, and Ba and ∼0.2 dex depletion of Cu in young clusters. Finally, we present a restricted three-process model where we add an asymptotic giant branch star (AGB) component to better fit Ba and Y. With the addition of the third process, we identify a population of stars, preferentially young, that have much higher AGB enrichment than expected from their SNIa enrichment. 
    more » « less
  6. Abstract We report the first direct measurement of the helium isotope ratio, 3 He/ 4 He, outside of the Local Interstellar Cloud, as part of science-verification observations with the upgraded CRyogenic InfraRed Echelle Spectrograph. Our determination of 3 He/ 4 He is based on metastable He i * absorption along the line of sight toward Θ 2 A Ori in the Orion Nebula. We measure a value 3 He/ 4 He = (1.77 ± 0.13) × 10 −4 , which is just ∼40% above the primordial relative abundance of these isotopes, assuming the Standard Model of particle physics and cosmology, ( 3 He/ 4 He) p = (1.257 ± 0.017) × 10 −4 . We calculate a suite of galactic chemical evolution simulations to study the Galactic build up of these isotopes, using the yields from Limongi & Chieffi for stars in the mass range M = 8–100 M ⊙ and Lagarde et al. for M = 0.8–8 M ⊙ . We find that these simulations simultaneously reproduce the Orion and protosolar 3 He/ 4 He values if the calculations are initialized with a primordial ratio 3 He / 4 He p = ( 1.043 ± 0.089 ) × 10 − 4 . Even though the quoted error does not include the model uncertainty, this determination agrees with the Standard Model value to within ∼2 σ . We also use the present-day Galactic abundance of deuterium (D/H), helium (He/H), and 3 He/ 4 He to infer an empirical limit on the primordial 3 He abundance, 3 He / H p ≤ ( 1.09 ± 0.18 ) × 10 − 5 , which also agrees with the Standard Model value. We point out that it is becoming increasingly difficult to explain the discrepant primordial 7 Li/H abundance with nonstandard physics, without breaking the remarkable simultaneous agreement of three primordial element ratios (D/H, 4 He/H, and 3 He/ 4 He) with the Standard Model values. 
    more » « less
  7. ABSTRACT We test the hypothesis that the observed first-peak (Sr, Y, Zr) and second-peak (Ba) s-process elemental abundances in low-metallicity Milky Way stars, and the abundances of the elements Mo and Ru, can be explained by a pervasive r-process contribution originating in neutrino-driven winds from highly magnetic and rapidly rotating proto-neutron stars (proto-NSs). We construct chemical evolution models that incorporate recent calculations of proto-NS yields in addition to contributions from asymptotic giant branch stars, Type Ia supernovae, and two alternative sets of yields for massive star winds and core-collapse supernovae. For non-rotating massive star yields from either set, models without proto-NS winds underpredict the observed s-process peak abundances by 0.3–$1\, \text{dex}$ at low metallicity, and they severely underpredict Mo and Ru at all metallicities. Models incorporating wind yields from proto-NSs with spin periods P ∼ 2–$5\, \text{ms}$ fit the observed trends for all these elements well. Alternatively, models omitting proto-NS winds but adopting yields of rapidly rotating massive stars, with vrot between 150 and $300\, \text{km}\, \text{s}^{-1}$, can explain the observed abundance levels reasonably well for [Fe/H] < −2. These models overpredict [Sr/Fe] and [Mo/Fe] at higher metallicities, but with a tuned dependence of vrot on stellar metallicity they might achieve an acceptable fit at all [Fe/H]. If many proto-NSs are born with strong magnetic fields and short spin periods, then their neutrino-driven winds provide a natural source for Sr, Y, Zr, Mo, Ru, and Ba in low-metallicity stellar populations. Conversely, spherical winds from unmagnetized proto-NSs overproduce the observed Sr, Y, and Zr abundances by a large factor. 
    more » « less
  8. ABSTRACT We develop a hybrid model of galactic chemical evolution that combines a multiring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradient in outflow mass loading tuned to reproduce the observed metallicity gradient. With this approach, the model reproduces many of the striking qualitative features of the Milky Way disc’s abundance structure: (i) the dependence of the [O/Fe]–[Fe/H] distribution on radius Rgal and mid-plane distance |z|; (ii) the changing shapes of the [O/H] and [Fe/H] distributions with Rgal and |z|; (iii) a broad distribution of [O/Fe] at sub-solar metallicity and changes in the [O/Fe] distribution with Rgal, |z|, and [Fe/H]; (iv) a tight correlation between [O/Fe] and stellar age for [O/Fe] > 0.1; (v) a population of young and intermediate-age α-enhanced stars caused by migration-induced variability in the Type Ia supernova rate; (vi) non-monotonic age–[O/H] and age–[Fe/H] relations, with large scatter and a median age of ∼4 Gyr near solar metallicity. Observationally motivated models with an enhanced star formation rate ∼2 Gyr ago improve agreement with the observed age–[Fe/H] and age–[O/H] relations, but worsen agreement with the observed age–[O/Fe] relation. None of our models predict an [O/Fe] distribution with the distinct bimodality seen in the observations, suggesting that more dramatic evolutionary pathways are required. All code and tables used for our models are publicly available through the Versatile Integrator for Chemical Evolution (VICE; https://pypi.org/project/vice). 
    more » « less
  9. null (Ed.)
    ABSTRACT We investigate the impact of bursts in star formation on the predictions of one-zone chemical evolution models, adopting oxygen (O), iron (Fe), and strontium (Sr), as representative α, iron-peak, and s-process elements, respectively. To this end, we develop and make use of the Versatile Integrator for Chemical Evolution (VICE), a python package designed to handle flexible user-specified evolutionary parameters. Starbursts driven by a temporary boost of gas accretion rate create loops in [O/Fe]–[Fe/H] evolutionary tracks and a peak in the stellar [O/Fe] distribution at intermediate values. Bursts driven by a temporary boost of star formation efficiency have similar effects, and they also produce a population of α-deficient stars during the depressed star formation phase following the burst. This α-deficient population is more prominent if the outflow rate is tied to a time-averaged star formation rate (SFR) instead of the instantaneous SFR. Theoretical models of Sr production predict a strong metallicity dependence of supernova and asymptotic giant branch star yields, though comparison to data suggests an additional, nearly metallicity-independent source. Evolution of [Sr/Fe] and [Sr/O] during a starburst is complex because of this metallicity dependence and the multiple time-scales at play. Moderate amplitude (10–20 per cent) sinusoidal oscillations in SFR produce loops in [O/Fe]–[Fe/H] tracks and multiple peaks in [O/Fe] distributions, a potential source of intrinsic scatter in observed sequences. We investigate the impact of a factor ∼2 enhancement of Galactic star formation ∼2 Gyr ago, as suggested by some recent observations. VICE is publicly available at <http://pypi.org/project/vice/>. 
    more » « less